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The pseudogap refers to an enigmatic state of matter with unusual physical properties found below a 
characteristic temperature T* in hole-doped high-temperature superconductors. Determining T* is critical 
for understanding this state. Here we study the simplest model of correlated electron systems, the Hubbard 
model, with cluster dynamical mean-field theory to find out whether the pseudogap can occur solely because 
of strong coupling physics and short nonlocal correlations. We find that the pseudogap characteristic 
temperature T* is a sharp crossover between different dynamical regimes along a line of thermodynamic 
anomalies that appears above a first-order phase transition, the Widom line. The Widom line emanating 
from the critical endpoint of a first-order transition is thus the organizing principle for the pseudogap phase 
diagram of the cuprates. No additional broken symmetry is necessary to explain the phenomenon. Broken 
symmetry states appear in the pseudogap and not the other way around. 
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The phase diagram of hole-doped high-temperature superconductors remains puzzling. A state of matter with 
unusual physical properties, dubbed "the pseudogap", is found below a characteristic temperature T* in a 
doped Mott insulator. Since the superconducting state is born out of the pseudogap over much of the phase 
diagram, the nature of the pseudogap is a fundamental issue in the field and it is under intense theoretical 1 and 
experimental 2 4 scrutiny. 

A pseudogap can occur because of disorder-broadened long-range ordered phases of Ising type 5 or because of 
fluctuating precursors to a long-range ordered phase that would appear only at T = 0 because of the Mermin- 
Wagner theorem 6 . The phenomenology of the pseudogap in the less strongly-coupled electron-doped cuprates 7,8 
differs from that of hole-doped cuprates and is consistent with fluctuating precursors 9 . The above two generic 
mechanisms have in common that they can occur even when interactions between electrons are not strong 
enough to lead to a Mott insulator at half-filling. Various broken symmetry states, such as spin-charge density 
wave, have been linked to the pseudogap in some cuprate families 2 4,9 but not in all of them. Yet the pseudogap is a 
generic feature of all hole-doped cuprates. It is also possible that broken symmetries are different phenomena, and 
in particular only a consequence, and not the origin, of the pseudogap, as suggested recently 10 . 

On the other hand, the Mott phenomenon, a blocking of charge transport because of strong electronic 
repulsion, is ubiquitous in hole-doped cuprates 11 and it is appropriate to ask whether there is a third generic 
mechanism for the pseudogap that is associated purely with Mott physics in two dimensions 7 . This mechanism 
appears only if the Coulomb repulsion is strong enough 7 to turn the system into a Mott insulator at zero doping, as 
observed in hole-doped cuprates. No broken translational or rotational symmetry is needed, although broken 
symmetry may also occur in certain cases 5 . In previous work 12,13 , we identified a first-order transition at finite 
doping between two different metals. Here we show that the characteristic temperature T* is an unexpected 
example of a phenomenon observed in fluids 14 , namely a sharp crossover between different dynamical regimes 
along a line of thermodynamic anomalies that appears above that first-order phase transition, the Widom line 14 . 

Results 

Widom line in doped Mott insulator. To investigate the formation of the pseudogap upon doping a Mott 
insulator, we study the competition between nearest neighbors hopping t and screened Coulomb repulsion U 
embodied in the two-dimensional Hubbard model. For the cuprates, the clear momentum dependence of the self- 
energy observed in photoemission (ARPES) forces one to use cluster extensions 15,16 of dynamical mean-field 
theory 17 . It is known that four sites 18 suffices to reproduce qualitatively the ARPES spectrum. Larger cluster sizes 19 
will improve momentum resolutions and change quantitative details but should not remove first-order 
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transitions, as has been verified at half-filling 20 . As in any cluster mean- 
field theory (dynamical or not), we can study the normal-state phase of 
the model by suppressing long-range magnetic, superconducting or 
other types of order while retaining the short-range correlations. 
Previous works based on this theoretical framework have shown that 
a pseudogap appears close to a Mott insulator, as a result of short-range 
correlations 718,21,22 . But it is only with recent theoretical and com- 
putational advances 23 that detailed results in the low temperature 
region of the phase diagram can be obtained. 

Recently, we used cellular dynamical mean- field theory on a 2X2 
plaquette to demonstrate that a first-order transition inhabits the 
finite doping part of the normal-state phase diagram of the model 12,13 . 
That first-order transition occurs along a coexistence line between 
two types of metals and ends at a critical point (T p ,fi p ) (see T — 
phase diagram in Fig. la). Here we show that the pseudogap is the 
low-doping phase whose characteristic temperature T* can be inter- 
preted as the Widom line. 

The Widom line is defined as the line where the maxima of dif- 
ferent thermodynamic response functions touch each other asymp- 
totically as one approaches the critical point 14 . The maxima become 
more pronounced on approaching T p , diverging at T p . Fig. lb shows 
the chemical potential dependence of the charge compressibility K = 
\ln 2 (dnldp) T for several temperatures above T p . Far above T p , k 
develops a maximum that increases and moves to higher doping with 
decreasing temperatures, indicating that the charge compressibility 
diverges at the critical point. The loci of the charge compressibility 
maxima, maxj^K, are shown in Fig. la and give an estimate of the 
Widom line. 

Crossing of the Widom line involves drastic changes in the 
dynamics of the system, as indicated in previous investigations on 
the phase diagram of fluids 14,24 . Similarly, here we show that the 
pseudogap, as seen in the single-particle density of states and in 
the zero frequency spin susceptibility, arises at high temperature 
from crossing the Widom line that radiates out of the critical point. 



Identification of T* on the basis of the local density of states. First 
we study the development of the pseudogap in the local density of 
states A(u>), which can be accessed by tunnelling or photoemission 
spectroscopy, along paths at constant temperature or at constant 
doping. Figs. 2(a,c) show the evolution of the density of states with 
doping at a fixed temperature above and below T p respectively. In the 
Mott insulating state, at 8 = 0, the density of states consists of lower 
and upper Hubbard bands separated by a correlation gap. Upon hole 
doping, there is a dramatic transfer of spectral weight from high to 
low frequency, as a consequence of strong electronic correlations. 
The low frequency part of A(io) develops a pseudogap, i.e. a 
depression in spectral weight, between a peak just below the Fermi 
level and a peak above the Fermi level. The pseudogap exhibits a two- 
peak profile that is highly asymmetric 22 , reflecting the large particle- 
hole asymmetry observed experimentally 25 . Upon increasing doping, 
the particle-hole asymmetry decreases, the spectral weight inside the 
pseudogap gradually fills in, and the distance between the two peaks 
slightly decreases. Finally, at a temperature-dependent doping, the 
pseudogap disappears and a rather broad peak appears in the density 
of states which narrows with increased doping. The change from 
pseudogap to correlated Fermi liquid behavior occurs either by a 
first-order transition when T < T p , with a discontinuous change in 
A(w), or by a crossover when T > T p . We obtain T*{5) from the 
inflection point in A(co = 0)(<5) at finite doping (see Figs. 2(b,d)). 

Figs. 2(e,g) show the temperature evolution of the density of states 
at constant doping. In Fig. 2(g) above 8 p (below fi p ), there is no 
pseudogap, only a peak around cu = 0 that broadens with increasing 
temperature 13 ' 22 . By contrast, in Fig. 2(e) below 8 p (above pi p ), the 
effect of increasing the temperature is to gradually fill the pseudogap, 
without decreasing the peak-to-peak distance. We identify the dis- 
appearance of the pseudogap in the spectrum by the inflection in A(co 
= 0)(T) (see Fig. 2(f) and Supplementary Fig. SI). We will discuss 
later the evolution with doping of X*, but the absence of pseudogap 
for S > 8 p , already allows us to conclude that T* appears as a bridge 
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Figure 1 | Phase diagram and Widom line, (a), temperature T - chemical potential /< phase diagram of the two-dimensional Hubbard model in the 
normal-state obtained by cellular dynamical mean-field theory. We take units where a = t = %= ks = 1. Here, U= 6.2 > Umit ~ 5.95 necessary to create 
a Mott insulator at half filling (;i = 0). The phase transition between two metals, the pseudogap and the correlated Fermi liquid, occurs along the 
coexistence line (dashed green line) that ends at the critical point (p. p , T p ) = ( — 0.707, 1/65) (orange circle). The coexistence line is bounded by the 
spinodals (blue lines with squares). The phase boundaries are obtained from scans at constant Tmonitoring the behavior of the occupation n versus /Z 2,13 . 
The onset of the Mott insulator (black line with squares) is defined by a plateau in the occupation atn= 1. The compressible, metallic, regions have n < 1. 
Below T p , the particle occupation n changes discontinuously at the spinodal lines, indicating the first-order transition. The green coexistence line and the 
extrapolations to T = 0 are a guide to the eye. The Widom line, i.e. the line where the maxima of different response functions converge, extends the 
coexistence line above T p . We estimate this line by the loci of charge compressibility maxima max,,K, where k = \l t?{dnl d\i) T (red line with triangles), (b) 
Semi-logarithmic plot of the charge compressibility k versus chemical potential ;i for several temperatures above Tp, obtained by numerical derivative of 
the filling with respect to the chemical potential. The loci of maxima of K are plotted in the T—fi plane in (a). The value at the maximum increases as T 
decreases, indicating the divergence of K at the critical point. This is analog to the diverging charge compressibility at the Mott critical point 33 or at the 
liquid-gas transition. 
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Figure 2 | T* from density of states. (a,c,e,g), A(co) at U = 6.2 obtained from maximum entropy analytic continuation of continuous-time quantum 
Monte Carlo data. (a,c), A(co) for different dopings at fixed temperature, (a) T = 1/25 > T p and (c ) T = 1/100 < T p respectively. In the latter case, 
associated with the first-order transition at finite doping there are hysteresis effects in the plane, where the transition appears as a coexistence interval 
between two metallic phases. Orange and magenta lines are A(oj) obtained in the coexistence region: they have the same chemical potential ;i= —0.725, 
but different filling and different low-energy spectrum, one with a pseudogap and the other one with a rather broad peak at oj = 0. (e,g) , density of states at 
fixed doping for various temperatures, (e) is for 5 = 0.02 < <5 p , where a pseudogap opens up with decreasing temperature, (g) is for <S = 0.08 > d p where a 
pseudogap is absent and a narrow peak develops as T decreases. (b,d,f,h), density of states at the Fermi level A(a> = 0) as a function of doping or 
temperature. Here, data are obtained from the extrapolated value of the imaginary part of the local cluster Green's function — l/7ilm G(co„ — » 0), which 
does not require analytical continuation. The inflection point of these curves as a function of /.i or T, indicated by a red circle, is our estimate of T*. On the 
upper horizontal axis we convert into physical units by using t = 0.35 eV. 



between the Mott insulator and the first-order transition at finite 
doping. 

Pseudogap versus Mott insulator. Next we discuss the relationship 
between the pseudogap and the Mott insulator. The pseudogap is 
linked to Mott physics because it opens up only above the critical 
17 for the Mott transition, in the metal near the Mott insulator. A 
pseudogap can occur below this threshold due to long-wavelength 
antiferromagnetic fluctuations, but that is different physics 6 occurring 
at a different energy scale. Here, only short-range spin correlations 
are involved, as observed experimentally 26 in YBa 2 Cu 4 0 8 at the 
pseudogap temperature. The first-order transition at finite doping, 
which is the terminus of the pseudogap phase, is linked as well to 
Mott physics because in the (17, T, (i) phase diagram, it emerges out of 
the Mott endpoint at half-filling, and progressively moves away from 
half-filling with increasing l/ 12,13 . 

The pseudogap inherits many features from the parent Mott insu- 
lator: in both phases the electrons are bound into short-range singlets 



because of the superexchange mechanism, reminiscent of the res- 
onating valence bond state 27 . Fig. 3a shows the temperature evolution 
of the probability obtained from the largest diagonal elements of the 
reduced density matrix on a 2 X 2 cluster 28 , which provide direct 
access to spin correlations (see also Supplementary Fig. S2). Below S p , 
upon decreasing temperature, the singlet is gaining weight at the 
expense of the triplet indicating a reduction of spin fluctuations. 
Because of singlet formation, the spin susceptibility y{T) drops below 
a characteristic temperature, as shown in Fig. 3b and as found in 
experiments 29 . The inflection point in y{T) defines a I* that moves 
to lower temperatures as the doping increases and approaches T„ 
as 3 — » S p . 

Despite the magnetic behavior similar to the Mott insulator phase, 
the pseudogap phase is a new state of matter. At T = 0 it appears to be 
separated from the Mott insulator by a second-order transition 13 . 
With decreasing temperature, neither the value of the density of 
states at the Fermi level (Fig. 2(f)), nor the spin susceptibility 
(Fig. 3b), extrapolate to zero. The T* extrapolated to i5 = 0 is not 
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Figure 3 | T* from spin susceptibility and plaquette eigenstates. (a), probability of the following plaquette eigenstates as a function of temperature, for 
several values of doping: the singlet N = 4, S = 0, K = (0, 0)) and the triplet N = 4, S = l,K = (n, n)) (squares and triangles respectively), where N, S, K 
are the number of electrons, the total spin and the cluster momentum of the plaquette eigenstate. (b), the zero-frequency spin susceptibility 
x{T) = J 0 (S z (t)S z (0))<?t as a function of temperature for several values of doping. S z is the projection of the total spin of the plaquette along the z 
direction. The inflection point of these curves as a function of T, marked by a solid symbol, is our estimate of I*. Data are for U = 6.2. 



related to the opening of the Mott gap (Fig. 4). Finally, the peak to 
peak distance for the pseudogap does not extrapolate, as S —* 0, to the 
Mott gap, just as is observed in experiment 30 . 




5 

Figure 4 | Pseudogap temperature T* along the Widom line. T-S phase 
diagram for U = 6.2 (see Supplementary Fig. S3 for the T — fi phase 
diagram). The pseudogap temperature scale T* is computed from the 
inflection point along paths at constant S of A(oj = 0)(line with green 
circles), x(T) (magenta circles) and Prob sing i < . t (T) (cyan circles). The 
T* line appears as soon as the Mott insulator is doped and joins the critical 
endpoint ( T p , d p ) (orange circle) of a first-order transition away from half- 
filling, closely following the Widom line (red line with triangles) that 
emanates from the critical endpoint. The critical point moves to large 
doping and low temperatures with increasing U 12,13 . The hatched region 
corresponds to the instability region bounded by the spinodals. 
Extrapolations to T = 0 are a guide to the eye. On the right vertical axis we 
convert into physical units by using t = 0.35 eV. 



T* as the Widom line. Now we move to the relationship between the 
pseudogap and the Widom line emanating from the critical point at 
finite doping. Fig. 4 shows in the T — S plane the doping evolution of 
the various T*, identified above as inflection points in A(a> = 0) 
along constant T or constant 8 paths, and as inflection points in 
X(T) and Prob sing i et (T) at constant S. The different T* lines for the 
pseudogap move closer to each other along the Widom line as we 
approach the critical endpoint. The interrelation between T* and 
Widom line is our main finding. Therefore our work shows that 
the dynamic crossover associated with the buildup of the pseu- 
dogap is concomitant with a crossover in the thermodynamic 
quantities, as observed in supercritical fluids. The organizing 
principle of these phenomena is the Widom line. One can thus 
interpret T* as the Widom line, or, equivalently, consider the 
Widom line as a thermodynamic signal of T*. Our results suggests 
that all indicators (both thermodynamic and dynamic) of the 
pseudogap temperature scale T* should approach each other with 
increasing doping, joining a critical endpoint of a first-order 
transition, which thereby appears as the source of anomalous 
behavior. In this view, pseudogap and strongly correlated Fermi 
liquid are separated from each other at low temperature by a first- 
order transition and are thus two distinct states of matter, just as 
liquid and gas are two distinct states, or phases. 

Discussion 

Common theories to explain the pseudogap phase include the pres- 
ence of rotational and/or spatial broken-symmetry phases as an 
essential ingredient. By contrast, in our approach the pseudogap is 
a consequence of large screened Coulomb repulsion that leads to 
strong singlet correlations in two dimensions reminiscent of resonat- 
ing valence bond physics 27 . Competing phases are not necessary to 
obtain a pseudogap. The pseudogap phase can however be unstable 
to such phases. We therefore provide a new and generic mechanism 
for the pseudogap in doped Mott insulators, according to which the 
pseudogap state is a new state of matter, whose characteristic tem- 
perature T* corresponds to the Widom line arising above a first- 
order transition. 

The Mott transition is often masked by broken-symmetry states. 
Similarly, our finite-doping transition is masked by the supercon- 
ducting phase for instance 31 . Nevertheless the rapid crossover be- 
tween pseudogap and metallic phases observed above the broken- 
symmetry states is accessible and can be controlled by the Widom 
line. Such rapid change in dynamics is a hallmark of the Widom 
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line 14,24 and it is consistent with the strong coupling nature and the 
observed phenomenology of the pseudogap in the vicinity of T* in 
hole-doped high-temperature superconductors. 

From a broader perspective, our work brings the conceptual 
framework of the Widom line, recently developed in the context of 
fluids 14,24 , to a completely different state of matter, the electronic 
fluid, suggesting its unexpected generality. We recall the strong 
impact that resulted from bringing in the field of electronic prop- 
erties of solids the well known concepts of smectic and nematic order 
developed earlier in the field of liquid crystals. It is tempting to argue 
that the same fate awaits the Widom line. 

Methods 

Our results are based on the cellular dynamical mean-field theory (CDMFT) 15 '' 6 
solution of the two-dimensional Hubbard model on the square lattice, 

H = - E ><>^ c i° +U H ("* - 0 (*j - i) - f E C 1 ) 

where c ia and operators destroy and create electrons on site i with spin <j, and 
n iG — c i(J . t is the hopping amplitude between nearest neighbors, U the energy cost 
of double occupation at each site and [X the chemical potential. CDMFT isolates a 
cluster of lattice sites, here a 2X2 plaquette, and replaces the missing lattice envir- 
onment by a bath of non- interacting electrons which is self-consistently determined. 
The cluster in a bath problem is solved by a continuous-time Quantum Monte Carlo 
summation of all diagrams obtained by expanding the partition function in powers of 
the hybridization between bath and cluster 23 ' 28 . The size of the plaquette is large 
enough to be consistent with the experimental observation 26 that at T*, in hole-doped 
cuprates, the antiferromagnetic correlation length is one or two lattice spacings. 

The value of Coulomb interaction U — 6.2 t is larger than the critical threshold 
t7 MIT ~ 5.95 t necessary to obtain a Mott insulator at half-filling 1213 and is chosen 
such that the pseudogap critical temperature T p (U) is accessible with our method. We 
carry out simulations at constant temperature for several values of fx and at constant 
doping 3 = 1 — n for several temperatures. Critical slowing down is a widespread and 
standard signal that the system is approaching a critical threshold 32 and appears in 
our simulations along the Widom line close to the critical point. To obtain reliable 
results, we perform up to 10 7 Monte Carlo sweeps, averaged over 64 processors, and 
hundreds of CDMFT iterations. The typical error on n is of order 10~ 5 . 
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